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Abstract 

Measurement of the distribution of hadronic energy in the final state in 
deeply inelastic electron scattering at HERA can provide a good test of our 
understanding of perturbative QCD. For this purpose, we consider the energy 
distribution function, which can be computed without needing final state 
parton fragmentation functions. We compute this distribution function for 
finite transverse momentum qj- at order a^, and use the results to sum the 
perturbation series to obtain a result valid for both large and small values of 
transverse momentum. 
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Typeset using REVT^ 



I. INTRODUCTION 



This paper concerns the energy distribution in the final state of deeply inelastic lepton 
scattering. Using a naive parton model, one would predict that the scattered parton appears 
as a single narrow jet at a certain angle (6'^,, 0^,) in the detector. Taking hard QCD interac- 
tions into account, one predicts a much richer structure for the final state energy distribution. 
In a previous paper (henceforth referred to as I), we investigated this structure, using an 
energy distribution function defined in analogy to the energy-energy correlation function in 
e~^e~ annihilation, i'i We studied this energy distribution as a function of angle {0b,4>b) in 
the detector in the region not too near to the direction {6^, 0*). In this region, simple QCD 
perturbation theory is applicable, and we presented calculations at order a^. In this paper, 
we extend the analysis to the region of {9b, (ps) near to . Here, multiple soft gluon 

radiation is important. Thus we use a summation of perturbation theory. 



A. The energy distribution function 

There is extensive literature on semi-inclusive deeply inelastic scattering; ffllJ a brief 
history and complete set of references can be found in paper I. We begin here with a concise 
review of how the energy distribution function is defined, and then discuss how we sum the 
contributions that are important in the region {6b,4'b) — {d*,4>*) to obtain a result which 
is valid for all values of {9b, (t>B)- 

The reaction that we study ise+A — > e + B + Xat the HERA electron-proton collider. 
Let us describe the particles by their energies and angles in the HERA laboratory frame, 
with the positive z-axis chosen in the direction of the proton beam and the negative z-axis 
in the direction of the electron beam. In completely inclusive deeply inelastic scattering, one 
measures only E' and 6'', the energy and angle of the scattered electron. In the semi-inclusive 
case studied in this paper, one also measures some basic features of the hadronic final state. 
In principle, one can measure the energy Eb and the angles {Ob, <Pb) of the outgoing hadron 
B. However, it is much simpler to perform a purely calorimetric measurement, in which 
only the total energy coming into a calorimeter cell at angles {9b,<Pb) is measured. This 
calorimetric measurement gives the energy distribution 



= y 1 dEB{l- cose b)Eb ^^(e + A^B + X) 

dE' dcosO' dcosOB dcpB b dcosO' dEB dcosOB dcpB 



The sum runs over all species of produced hadrons B. We have included a factor (1 — cosOb) 
in the definition because this factor is part of the Lorentz invariant dot product Pa,^iPb = 
EaEb{1- cos Ob). 

Notice that dH measures the distribution of energy in the final state as a function of 
angle without asking how that energy is split into individual hadrons moving in the same 
direction. El For this reason, the theoretical expression for dT, will not involve parton decay 
functions that describe how partons decay into hadrons. 
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B. Partonic variables and their relation to HERA lab frame 

At the Born level, the hard scattering process for the reaction is electron + quark — »• 
electron + quark by means of virtual photon or Zq exchange. At order a^, one can have 
virtual corrections to the Born graph. In addition, one can have processes in which there are 
two scattered partons in the final state. Then the initial parton can be either a quark (or 
antiquark) or a gluon, while the observed hadron can come from the decay of either of the 
final state partons. Some of these possibihties are illustrated in Fig. |I|, Fig. ||, and Fig. |[ 

Let us consider the effect of the emission of the additional, unobserved, "bremsstrahlung" 
parton. We can define the part of the vector boson momentum g'^ that is transverse to the 
momentum of the incoming hadron and to the momentum of the outgoing hadron P^. 
One merely subtracts from q^ its projections along P^ and Pg (taking P^ = P§ = 0), 
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FIG. 3. Feynman diagrams for gluon initiated process with a quark jet observed. The observed 
parton is the upper hne, indicated with a dot. 
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(2) 



We let Qt = [—Qt ' Qt^] represent the magnitude of the transverse momentum. It is qt 
that is analogous to the transverse momentum of produced W's and Z's or lepton pairs in 
the Drell Yan process. In the naive parton model, there are no bremsstrahlung partons and 
all parton momenta are exactly collinear with the corresponding hadron momenta, so one 
has Qt = 0. At order unobserved parton emission allows qt to be nonzero. 

In order to properly describe the parton kinematics we need four more variables besides 
Qt- Two are the standard variables for deeply inelastic scattering, = —q'^q^ and x = 
I (2g ■ Pa)- The third is a momentum fraction for the outgoing hadron 5, 



Pb-Pa 



q-PA 



(3) 



(Thus the integration over the energy of hadron B in definition (|l|) of the energy distribution 
is equivalent to an integration over 2;.) The fourth variable is an azimuthal angle 0. To define 
0, we choose a frame, called the hadron frame. Fig. ^ in which the incoming hadron A has 
its three-momentum along the positive 2;-axis and the virtual photon four-momentum 
q^ lies along the negative 2;-axis. Then as long as 7^ 0, hadron B has some transverse 
momentum, and we align the x- and y-axes so that PJ > and Pg = 0. We now define 
as the azimuthal angle of the incoming lepton in the hadron frame. These variables are 
described more fully in paper I, and relevant formulas are given in the Appendix of this 
paper. 

The variables and can be translated to the observables of the HERA lab frame. 
Fig. |[ In the naive parton model, the outgoing hadron B (along with all the other hadrons 
arising from the decay of the struck quark) emerges in the plane defined by the incoming 
and outgoing electrons at a precisely defined angle {9^,(f)^), which can be computed from 
the incoming particle momenta and the momentum of the scattered electron. The point 
qT = corresponds to {Ob,4>b) = (6'*,0*)- We choose our x-axis such that 0^, = 0. Lines of 
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Pb +x 




FIG. 4. The hadron frame. The initial hadron Pa hes along the positive z-axis, and the vector 
boson q lies along the negative z-axis. The next-to-leading order QCD corrections can give the 
final state hadron Pb a non-zero x-component. 




X 



FIG. 5. The HERA lab frame for the process: e^{l) + A{Pa) e-{£') + B{Pb) + X. The 
final state leption e~{i') lies in the x-z-plane, and the final state hadron B{Pb) has a non-zero 
2/-component (j)B is non-zero. 

constant positive qx are curves in the {9b, (ps) plane that encircle the point (6',,,0^,). Lines 
of constant (j) radiate out of the point (6'*,0*), crossing the lines of constant qt- This is 
illustrated in Fig. ||. The precise formulas for the map relating {Ob,4>b) and {qT,(p) are 
given in the Appendix. 

In paper I and in this paper, we find it convenient to convert from the laboratory frame 
variables {E',9'} of the scattered lepton and {9b, (ps} of the observed hadron to {x,Q^} for 
the lepton and {qt, 0} for the observed hadron. We also convert from Eb to z. With this 
change of variables, Eq. (0) becomes 

= y[dzz(^] . (4) 

dx dQ"^ dqj' d4> J \2xEa) dx dQ^ dqj^ dcj) dz 



C. The Sudakov summation of logarithms of qt 

The main object of study in this paper is the distribution of energy as a function of qt 
for -C Q^. In paper I, we applied straightforward perturbation theory to analyze the 
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Q= 30 GeV x=0.02 



+71 



B 



-71 




■1 



Cos(Q^) 



■1 



FIG. 6. Contours in (pB and cos{9b) for Q = 30 GeV and x = 0.02. The circular rings are 
contours of constant qt in steps of 3 GeV, and the radial arcs are contours of constant (p in steps 
of tt/8. 



energy distribution in the region ~ Q^, and asiQ"^) <^ 1. Here there is a rich structure 
as a function of the angles that relate the hadron momenta to the lepton momenta. In fact, 
a complete description requires nine structure functions. 

When one examines the region Q^, one finds that the angular structure simplifies 

greatly. However, the dependence on g-p becomes richer than the dependence on q^p of the 
lowest order graphs. By summing the most important parts of graphs at arbitrarily high 
order, one finds a structure that is sensitive to the fact that QCD is a gauge theory. 

Briefly, the physical picture is as follows. At the Born level of deeply inelastic scat- 
tering, a quark in the incoming proton enters the scattering with momentum ^Pj^ that is 
precisely along the beam axis. This quark is scattered by a virtual photon, Z- or IV-boson. 
Its momentum ^P^ + q^^ is in a direction [6^, cf)^) that can be reconstructed by knowing the 
lepton momenta. However at higher orders of perturbation theory, the momentum of the 
final state parton is 



{ki -\- k2 -\- ■ ■ ■ -\- J 



(5) 



where the fcf are momenta of gluons that emitted in the process. In a renormalizable field 
theory, it is very easy to emit gluons that are nearly coUinear to either the initial or final 
parton directions. In addition, in a gauge theory such as QCD, it is very easy to emit gluons 
that are soft {k^ <^ Q). Each gluon emission displaces qx by a small amount, so that one 
may think of the parton as undergoing a random walk in the space of transverse momenta. 
With one gluon emission, one finds a cross section that is singular as 0: 



da a + blog{q^/Q^) 

oc as — 



dqj^ 



Or 



(6) 



At order the l/q^ singularity is multiplied by a polynomial in \og{q'^/Q'^) of order 2A^ — 1. 
This series sums to a function of that is peaked at g^ = but is not singular there. The 
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width of this distribution is much bigger than the 300 MeV that one would guess based on 
experience with soft hadronic physics. On the other hand the width is quite small compared 
to the hard momentum scale Q. 

Essentially this same physics has been studied in the two crossed versions of the process 
e + A —>■ e + B + X that can be studied at HERA. In electron-positron annihilation, 
e + e— > A + B + X , one looks at the energy-energy correlation function for hadrons A 
and B nearly back-to-back. iUll In A + B i + I + X, one studies the distribution 
of the lenton pair as a function of its transverse momentum qt with respect to the beam 
axis. 0O"ii The same analysis applies also to the distribution of the transverse momentum 
otWoiZ bosons produced in hadron colliders. 

From these studies, the following picture emerges. First, the leading logs (n = 2N — 1) 
can be summed to all orders, and dominate the perturbation theory in the region as(Q^) -C 1 
and as(Q^) ln^(Q^/g|.)^l. Unfortunately, most of the interesting physics, and most of the 
data, lie outside this region of validity of the leading logarithm approximation. Fortunately, 
one can go beyond the leading logarithm summation to obtain a result that is valid even 
when a,(Q2) In^Q^q^) is large. 011 

The plan for the remainder of the paper is as follows. In Sec. 2, we use our as calculation 
in paper I to calculate the asymptotic form of the energy distribution functions in the qt —>■ 
limit. In Sec. 3, we introduce the Sudakov form factor which sums the soft gluon radiation 
in the limit qx — > 0. In Sec. 4, we compare the asymptotic form of the energy distribution 
functions to extract the order as contributions to the perturbative coefficients A, B, and 
Sec. 5, we address the issue of matching the small qt region to the large qt region. 
In Sec. 6 we investigate the form of the non-perturbative corrections in the small qt region, 
and relate these to the Drell-Yan and e~^e~ processes. In Sec. 7, we review the principle 
steps in the calculation. In Sec. 8, we present results for the energy distribution functions 
throughout the full qt range. Conclusions are presented in Sec. 9, and the Appendix contains 
a set of relevant formulas. 



II. THE ENERGY DISTRIBUTION FUNCTIONS 

In this section we review the order as perturbative results of paper I in order to extract 
the terms in dl2 / dx dQ"^ dq^ dcj) that behave like l/q^ times logs as qt — > 0. In Sec. 3, we 
display the structure of / dx dQ^ dq^ dcf) with the Sudakov summation of logarithms. Then 
in Sec. 4, by comparing the summed form with the order as form of dTj/dx dQ"^ dq^d(j), we 
will be able to extract the coefficients that appear in the summed form. 



A. Energy Distribution Formulas 



The process we consider ise + A ^ e + B + X, and the fundamental formula for the 
energy distribution is: 

dH ^ 

^ ,Q2d2^^ = E E E ^o{Q';VuV2,j,j',k) rk{x,Q\q'T;j,j') . (7) 

^ "iT k=l Vi,V2 
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The hyperbolic boost angle, ip, that connects the natural hadron and lepton frame is given 



by [§ 



2x5 

cosh^ = — -1 , (8) 



and (f) is the azimuthal angle in the hadron frame. The nine angular functions Ak{ip, 4>) arise 
from hyperbolic D^[ijj,(l)) rotation matrices. The complete set of Ak{ip,4>) are listed in the 
Appendix, but the two we shall focus on are 

Ai{ilj,(f)) = 1 + cosh^(^) 

A(^,0) = 2 cosh(V^) . (9) 

We sum over the intermediate vector bosons {Vi,V2} = {7, Z"} or {ly^}, as appropriate, 
and we also sum over the initial and final partons, {j,j'}- The factor So(Q^; Vi, V2, j, j', k) 
contains the leptonic and partonic couplings, the boson propagators, and numerical factors; 
it is defined in the Appendix, Eq. (P^). It is the hadronic energy distribution functions, 
Tk{x,Q'^,q^;j,j'), that we shall calculate. 

If we expand the Tk in the form of perturbative coefficients convoluted with parton distri- 
bution functions, then two of the functions Tk, namely Fi and Fg, behave like log^ {q^/Q'^)/q^ 
with n > for qt — > 0. The others behave like or 1 times possible logarithms. In this 
paper we are interested in small qt behavior, so we concentrate our attention on Fi and Fg. 

What of the less singular structure functions F2, F3, F4, F5, F7, Fg and Fg? Fixed 
order perturbation theory is not applicable for the calculation of these F^ for small qt- 
We note, on the grounds of analyticity, that these Fj^ must be finite or, for certain k, 
vanish as qx —>■ 0, even though they have weak singularities in finite order perturbation 
theory. Our perturbative results in the region of moderate qt indicate that the fraction of 

/ dx dQ'^ dq^ dcp contributed by these F^ is small and dropping as qt decreases. We thus 
conclude that these contributions would be hard to detect experimentally for small qt- For 
this reason, we do not address the problem of summing perturbation theory for F2, F3, F4, 
F5, F7, Fs and Fg. 

Applying the methods of Refs. p8| , p!4| to deeply inelastic scattering, we write Fi in the 
form 

^l(x,Q^g^;J;/) = Ff-*(x,Q^g^;J■/)-Ff^"^(x,Q^g2.^.y)^^(^^g2^^2.^.^ _ 

(10) 

Here W{x,Q'^,q^;j,j') sums the singular terms to all orders, and contains the leading 
behavior of Fi as gr — > 0. Tf^^^{x,Q'^,q^;j,j') is simply Ti{x,Q'^,q^;j,j') evaluated at 
a finite order (aj for our purpose) in perturbation theory. Tf'^^"^{x,Q'^,q^;j,j') equals 
W{x,Q^,q^;j,j') truncated at a finite order of in perturbation theory. Specifically, if 
we expand W{x, Q^, q^]],]') in the form of perturbative coefficients convoluted with parton 
distribution functions, then the coefficients have the form of log" {q^ / Q'^) / q^ with n > 0. 
There are, by definition, no terms that behave like {q^/Q'^Y times possible logarithms for 
p > —1. Such terms exist in Fi, but they are associated with (rf*^*"* — Yf'^^"^) in Eq. (p!0[). 

The angular function Ai{ip,(l)) = 1 + cosh^ (■?/') that multiplies W in the small qx limit 
arises from the numerator factor 



8 



Drell-Yan DIS e+e" 

FIG. 7. Dominant integration regions leading to the non-perturbative contributions to (a) 
Drell-Yan, (b) DIS, and (c) e^e^. These three processes are related via a crossing symmetry. Ja 
and Jb represent the jet subgraphs associated with the collinear partons from hadron A{Pa) and 
B{Pb), respectively. S represents the subgraph of soft gluons and quarks which are connected to 
the rest of the process by soft gluons (but not soft quarks). The double-logarithms arise from Ja 
and Jb- 

TT{h/i.} TiitiY^BY} . (11) 

Here Tr{/7^/ 7,^} is associated with the lepton scattering, and the factor ■ ■ ■ ^ gives the 
Dirac structure of the hadronic part of the cut diagram in Fig. |^(b) in the limit — > 0. We 
will discuss Fig. |^ further in Sec. 0. 

The weak currents also contain 757^ terms. This gives the possibility of another angular 
function in the small qt limit. With the same limiting hadronic structure, I/'a ■ ■ ■ I/'b can 
have 

Tr{/757/7.} "TiiAl^r^BY], (12) 

which is proportional to the angular function AQi^^cf)) = 2cosh('?/') at qt = 0. (Note that 
both Ai{ip,(j)) and AQ{ip,(j)) are independent of the azimuthal angle 0.) Thus Tg has the 
structure 

^6(x,Q^g^;J,/) = *(x,Q^g|;J,/) -^f^"^(x,Q^g^;J,/) + (-l)W^(x, Q^ j,/), 

(13) 

with the same function^ W as in Eq. (p^. Again W contains the terms that behave like 
log"'(g|./Q^)/q'T perturbation theory, while (T^^''* — T^'^^™') contain the less singular terms. 
Our object now will be to study the small qx function W. 



^ The minus sign in front of W{x,Q'^,q^;j,j') in Eq. ( [T^ ) arises from our convention for the 
functions ^^(V','/') and couplings So((5^; Vi, V2, j, j', fc) that multiply Ti and Tg. 
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B. Parton Level Distributions 



The above hadronic process takes place via the partonic sub-process V{q) + a{ka) 
b{ki,) + X where V is an intermediate vector boson, and a and b denote parton species. The 
hadron structure function W{x, Q^, is related to a perturbatively calculable parton 

level structure function Wa{x,Q^,q^;j,f) via 

W{x, Q^, q^; j, j') = fa/A ®Wa^ "T X] fa/A{.i, 1^') Wa{x, , q^] j, f) , (14) 

with = / Pa ^^'^ ^ = ^/^a- Here fa/A is the MS parton distribution function. Note 
that the decay distribution function dB/b{ib,l^) is absent since we have used the extra J zdz 
and the sum over hadrons from the definition of the energy distribution to integrate out the 
ds/biCb, A*) via the momentum sum rule, 

E / dCUbdB/b{Cb,fi)^l ■ (15) 

B •' 

The partonic structure function, Wa{x,Q'^,q^;j,j'), is obtained by first computing the 
partonic tensor 

w''''{ka,kb,q)^^ E / d^xe-"^-^{ka,s\f{0)\kb,s';X){kb,s';X\f{0)\ka,s), (16) 



which is a matrix element of current operators. We then project out the appropriate angular 
component (c/., paper I), and extract the leading term in the qr ^ limit. Explicit 
calculation will show that these limits (up to overall factors) are identical for the projection 
of the 1 and 6 tensors. In the small qx limit, the energy distribution function is then given 
by: 

——^-^—c:,A^i^,4>) E J2^o{Q';V,,V2,j,j',i)J2fa/A{^,i^)®M^,Q',4;jJ') 

ax aqrp acp y^y^ ^j, „ 

-Ae{^,<f>) E J2^o{Q';V,,V2,j,j',6) ^ fa/A{^,f^)^Wa{x,Q',q^;j,f) 

Vi,V2 a 

+ plus terms less singular than 1/q^ . (17) 
Again, the relative minus sign is simply due to the definition of Aki'ip, (f>) and 

So(g';Vi,^,j,/,^). 



C. The Asymptotic Energy Distribution Functions 

We observe (from the results of paper I) that the perturbative Tf'^'^*{x,Q'^,q^;j,j') and 
r|'^''*(a;, Q^, Q'r; J, i') diverge as l/qx for qr — > 0. To identify the singular terms, we can 
expand the on-shell delta function for small qr using 
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27r 6 



MN2 
b ) 



2nx 



{■"(f)- 



x) 5(1 -z) + 



S{l-z) S(l-x) 



:i-x) 



'l-z] 



where the "+" -prescriptions is defined as usual by: 



(18) 



(19) 



Taking the qx ^ hmit for the results of paper I, we find the partonic energy distribution 
to be 



wt'y''\x,Q\q'^;j,f) 



167r^ as 



Sa,j 5(1-2) Cf 



2 In 



.4. 



l + x^ 

1 — X 



x' + il 



X] 



(20) 



where we use 5aj and 5a,g for the quark and gluon contributions, respectively. For conve- 
nience, we denote the asymptotic limit g-r of Wa by w^^^"^. 

In this limit, we can greatly simplify this expression by identifying the QCD splitting 
functions. We present the result for the hadronic structure function convoluted with the 
parton distributions, (c/., Eq. (|T^) ): 



^^^^"(x,Q^g|; 



167r^ as 
ql 



{fj/A{x) Ct 



+ fjlA ® Pq/q + fg/A ® Pq/g] 



(21) 



where ® represents a convolution in x. In the simple form above, it is easy to identify 
the separate contributions. The last two terms arise from the collinear singularities, and are 
proportional to the appropriate first order splitting kernel, Pq/q and Pq/g- It is the remaining 
term in which we are interested as these arise from the soft gluon processes. We note that 
pAsj/m jg (jggned such that the combination r^^''* — r^*^"^ has only logarithmic singularities 
as qT — * 0. 



III. SUDAKOV FORM FACTOR 

In this section, we display the structure of dTj / dx dQ"^ dq^ dcf) with the Sudakov summa- 
tion of logarithms. This provides the basis a formula that includes nonperturbative effects, 
developed in Sec. 6. In addition, in Sec. 4 we compare the summed form of this section with 
the order as form of dT, / dx dQ"^ dq^ d(j) from Sec. 2, in order to extract the coefficients that 
appear in the summed form. 
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A. Bessel Transform of Wa{x,Q'^,q^;j,j') 

It proves convenient to introduce a Fourier transform between transverse momentum 
space (qt) and impact parameter space (6), 

= / ^hUhqT)wa{x,Q\h''-j,j') , (22) 

as Wa{x,Q'^,b'^]j,j') will have a simple structure. Effectively, we make use of the renor- 
malization group equation to sum the logs of Q^, and gauge invariance to sum the logs of 
~ 1 /b. The Fourier transform also maps the qt singularities at the origin to the large b be- 
havior of Wa{x, Q^, b'^]j, j')] we will take advantage of this when we consider non-perturbative 
contributions. 



B. Sudakov Form Factor 



The structure function in impact parameter space Wa{x,Q^ ^b"^; j, j') has the factorized 
form: 



Wa{x,Q\b^-,j,j') = C'^^{xM J2 J dzzC^,fiz,bfi) e 



■5(6) 



(23) 



This form is from references and [14] applied to the DIS process, and generalized to 
include vector bosons other than the photon. The last exponential factor is the Sudakov 
form factor: 



S{b) 







|ln 













(24) 



The logarithm in the exponential is characteristic of the gauge theory. It arises from the 
soft gluon summation in QCD at the low transverse momentum -C Q^. The arbitrary 
constants {6*1,6*2} reflect the freedom in the choice of renormalization scale. We choose 
{C*i,C*2} to be 



C2 = 1 



(25) 
(26) 



The functions A, B and the hard scattering functions C*'s are simple power series in the 
strong coupling constant Us with numerical coefficients:^ 



^ Collins and Soper (CS) expand in powers of ag/iT, and Davies, Webber, and Stirling 
(DWS) expand in powers of as/(27r). We carry the extra factor of (2) explicitly to facilitate 
comparison between these references. 
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^(«^(/^))=E {^] An (27) 



N 



N=l 




N 



(28) 



oo 



N 



C^,^{xM = 6{l-x)6,^+Y. C]T\xM) {^[ (29) 



N=l 
oo 



N 



N=l 



The normalization has been chosen such that each hard scattering function C equals a 
(5-function at leading order. 



As noted in reference |TJ], in the limit Q ^ oo, all logarithms may be counted as being 
equally large. Therefore, to evaluate the cross section at — to an approximation of 
"degree A^," one must evaluate A to order a^^"*"^, B to order as^~^^, C^^ and C^^'^ to order 
and the j3 function order Us^^"^ ■ In particular, an extra order in A is necessary due to 



a 



N 



the extra logarithmic factor in Eq. (^J). For the present calculation, we evaluate A to order 
a^^, B to order a^^, C^^ and C^^'^ to order a^^, and the j3 function to order a^^. This 
yields the cross section to order for large g^, to order aj^ for small gr, and the cross 
section integrated over qx to 



C. Perturbative Expansion of the Sudakov Form Factor 

We can extract the Ai and Bi coefficients of the Sudakov factor by expanding 
'Wa{x,Q^,b'^;j,j') of Eq. ( P5D in a^, and comparing with the perturbative calculation of 
paper I. Here, we take a fixed momentum scale /xq in ctsifJ^o) as the running of as{fi) con- 
tributes only to higher orders. We can now compute the integral over //^ analytically to 
obtain: 

sin)) + B{as{^l)) 

(31) 

where 

(32) 

We expand the Sudakov exponential out to order a^^, 

e-^W ^ 1 - S{h) + , (33) 



S{h) 





In 




!ci/b^ /i^ 







A{c 



(2)7r 



^2 



In 
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and perform the Bessel transform of Wa{x, Q^, b'^;j,j') {cf., Eq. (p3D) to obtain the partonic 
structure function in momentum space: 



Jo Zir 



6aj6{l -x) + 



X 



l-S{b)+0{as 



(34) 



where we have used the first order expressions for C^^ (x, and Cj 
Finally, we integrate to obtain the 0{as^) terms for finite g^: 

167r^ a 



yOUT (N) 



z, fib). 



2Bi 



+ Sa,j Pqlq{x) + Pq/g{x) 

+ terms p..opo..tio„aUo^(,|)' 



(35) 



Here, we have use the fact that the renormalization group equation tells us the form of 
C^^^^\x, fib) and C'^^^^^\z, fib) must be a splitting kernel times log[/i6], plus a function 
independent of fi and b. Equivalently, for the hadronic structure function, we find: 



W{x,Q',q^;j,j' 



(It 



fj/A{x) 



^|lnm-21nr^ 
(2) 1 [q'^J \ 2C2 



2Bi 



+ fj/A ® Pq/q + fg/A ® Pq/g 

^ terms, „n.HoH4) 



(36) 



We will compare the first-order expansions in Eq. ( |35D and Eq. (|36| ) with the asymptotic 
limit of the perturbative calculations of Sec. |I| to extract the desired Ai and Bi coefficients. 

IV. COMPARING ASYMPTOTIC AND SUDAKOV CONTRIBUTIONS 

In this section, we compare the summed form of dT, / dx dQ"^ dq^ dip with the order 
form, and thus extract the coefficients that appear in the summed form. 

A. Extraction of A and B 



Comparing the expansion of the Sudakov expression [Eq. (pG])] with the asymptotic 
results [Eq. (PT|)], we obtain the order a], coefficients Ai and Si, 
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Ai = (2) Cf 
Bi = (2) 2 Cf In 



2C2 



.7i5-(3/4) 



With our particular choice of the arbitrary constants {Ci, C2} in Eq. 

A, = (2) Cf 
Bi = (2) 



I), we have: 



(37) 
(38) 

(39) 
(40) 



We find that the results for Ai and Bi obtained above are identical to those found in refer- 



ence 



^ for Drell-Yan production, as well as those found in reference [|T5| for e"*" e annihila- 
tion. This apparent crossing symmetry has been demonstrated at order by Trentadue. il 



In light of this result, we shall make use of the A2 coefficient & 

.2 in o . Ci 



(41) 



The extra order in the Ai expansion will compensate extra logarithm L which is not present 
for the Bi terms. 



B. Expansion of C^^ and C"^^^ 



C^^ and (7*^^^ terms are obtained by comparing the terms in the perturbative expansion 
proportional to SIqt) with the expanded summed form. Since the virtual graphs yield 
contributions only proportional to S^qr), they will only enter and C^^^ . The real 
graphs yield both zero and finite qr terms; therefore, they will contribute to both Aj, Sj, 
and the C^^ and C**^^^ coefficients. The calculation of the virtual graphs has been performed 
by Meng, and we make use of those results. 

We have defined the C{x,fib) coefficients such that at leading order, they are 

C'^r^'\^^.l^h) = 5,,5{l-z) 



C]^^'\x,,b)=C':r^'\z,,b) = ^ 



(42) 



(Here, j and k denote quarks and anti-quarks, and g denotes gluons.) At next to leading 
order, we find that C^^ '^^-^ (x, jjb) match those calculated by CS for the Drell-Yan process: o 



C]^^'\x, f,b) = 5,,|^(1 -x)+ P,/,{x) In (^A^ 
+ S{l-x] 



-Cf In' 



IT 



a A 



23 
12 



(43) 
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C;^^'\x,f^b) = -x) + P,/,(x) In (-^] . (44) 



The C'^'-^^ ^^\z , fxb) are simply those for e~^e as given in reference [|T5 



Cfr^'\^,f^b) = sJ^il -z) + P,i,{z) In 



+ 5(1 - z) 



-Cf In 



2 




(45) 



C,7^(^)(z,/x6) = f ?+P^/,(2) In ( 4 ) , (46) 



where we define A = 2e '^^ to simphfy the notation. Note that C^^ and C^^^ are only a 
function of the ratio C1/C2. 

C. Complete Expression 

Now that we have obtained Ai, A2, and Bi, we can substitute into equation Eq. (|24D 
to obtain the complete Sudakov contribution (including the full a<((/i) dependence). We 
choose to perform the fi integral analytically, as the Bessel transform would be prohibitively 
CPU intensive if we did not. To facilitate this computation we provide an integral table in 
Appendix |^ including all the necessary terms We are now ready to combine the separate 
parts of the calculation. 

V. MATCHING 

We now have computed the contributions to the energy distribution functions for the 
perturbative Ff'^'^* in paper I [Eq. ([T^ ], the summed (or Sudakov) W{x,Q'^,q^;j,j') in 
Eq. (^), and the asymptotic p^''^'" in Eq. (|2T|). We can simply assemble these pieces to 
form the total structure functions via: 

^l(x,Q^g^;J,/) = ^f-*(x,Q^g^;J,/)+ l^(x, g^ j,/) - rf ^"(x, Q^ g^; j,/) 
T,{x,Q\q'^-j,j') = T^^^\x,Q\q'r;j,j') + hl)W{x,Q\q'^-j,f)-Tt'^^^ 

(47) 

Here, F^'^''* and p^^^"^ are evaluated at order a^^, while W contains a summation of pertur- 
bation theory. In the limit ~^ 0, F^*^^'* and F^*^™ will cancel each other leaving W as we 
desire. In the limit qx — Q, W and F^*^™ will cancel to leading order in a^; however, the 
finite difference may not be negligible. To ensure that we recover the proper result (Ff ^''*) 
for large gr, we define the total energy distribution function (F^) to be: 



16 



0.2 0.4 0.6 0.8 1 

FIG. 8. The matching function T{qT/Q) vs. {qr/Q) for p = {1,3,5,7,9}. p = 1 is the top 
curve, and p = 9 is the bottom curve. 

^6(x,g^g^;J,/) = rf ^*(x,Q^g^._^. 

+ r (^l) { (-1) g^ g^; J, f) - rf ^"(x, g^ J, j')} , 

(48) 

where we introduce the arbitrary function 



^1^1=^^ ■ (49) 



Q 



The transition function T{q'j'/Q) serves to switch smoothly from the matched formulas to 
the perturbative formula, and p is an arbitrary parameter which determines the details of 
the matching. Fig. |] displays Ti^qrlQ) for a range of p values. We will choose p = 5 which 
ensures that ~ r^'e''* for g7^/Q,^0.4, a conservative value. 

VI. NON-PERTURBATIVE CONTRIBUTIONS 

In analogy with Eq. (|^) and Eq. (p3D , the Bessel transform of the hadronic structure 
function is defined as: 



W{x,Q\q'^;j,f) = / ^ e^'^' W{x,Q\b';j,f) . (50) 
When b is small, we have: 
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Wix,Q',b';j,j') = /'f E fa/A^^f^) Ci:!{xM E / dzzC^,f{zM e-'^'l (51) 

? a a' 

The perturbative calculation of W{x,Q'^,b'^;j,j') is not reliable for 6^1/A. However, the 
integration over b in Eq. (^) runs to infinitely large 6^, and the region 6^1/A is important 
for values of and of practical interest. In order to deal with the large 6^ region, we 
follow the method introduced in Refs. [HjQ- We define a value 6max such that we can 
consider perturbation theory to be reliable for b < femax- (In our numerical examples, we 
take l/6max = 2GeV.) Then we define a function 6* of b such that b^, ~ b for small b and 
b* < &max for all b: 



6* = , . (52) 

^1 + &V&Lx 



We define a version of W{x, Q^, b'^',j,j') for which perturbation theory is always reliable 
by W{x,Q^,bl;j,j'). Note that for small b, the difference between 14^(6*) and W{b) is 
negligible because b^ ~ b. Conversely, perturbation theory is always applicable for the 
calculation of W^(&*) because 6^, is small even when b is large. 

Next, we define a nonperturbative function exp(— S'np(&)) as the ratio of W{b) and Vr(6*): 

W{x,Q',b';j,f) = W{x,Q',bl-j,f) e-^-(-'«^^^^^-.^') . (53) 

Ultimately, we will have to use nonperturbative information to determine S^p{b). However, 
some important information is available to us. From Eq. (|5T|), we see that 

d\og[W{x,Q\b';j,f)] 

91ogQ2 ^ ^ 

is independent of x,j,j' and Q^. This result is derived in perturbation theory, but at 
arbitrary order, so we presume that it holds even beyond perturbation theory. Then 

^5Np(a:,g^&^;J,/) 

9b^Q^ ^^^^ 

is also independent of x,j,j', k and Q^. That is, S'np has the form 

SMx,Q',b';j,f) = \og(QyQl)gi{b) + ASM^,b'-j,j') . (56) 

(Here Qq is an arbitrary constant with dimensions of mass, inserted to keep the argument 
of the logarithm dimensionless.) Furthermore, in W, the x and j dependence occurs in a 
separate factor from the j' dependence. Thus the second term in Eq. ( |5^ ) above can be 
simplified to 



SM^,Q^b';j,f) = \og[QyQl)g,ib)+gAix,b';j)+gBib';f) . (57) 
(Recall, we have integrated over z.) 
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FIG. 9. Comparison of the non-perturbative function e 



,{b) 



vs. 



b for the Drell-Yan 



(Davies, Webber, and Stirhng ||l^) [upper line], and e'^e (Colhns and Soper 
for Q = 30GeV [Fig. (a)], and Q = 100 GeV [Fig. (b)]. 



[lower line]. 



Perturbation theory is not applicable for the calculation of the functions gi{b), gA{x, b'^', j) 
and gB{b^]j') for large b. For small b, perturbation theory tells us only that these functions 
approach as 6 — 0. This follows from Eq. (|5B[), and the fact that b^/b when 6^0. 
(See Ref. [Jl^ for further discussion.) Since we learn little from perturbation theory, we turn 
to non-perturbative sources of information. Fortunately, the analogous functions in e^e~ 
annihilation and in the Drell-Yan process have been fit using experimental results. [0011 

We therefore ask whether the functions gi{b), gA{x, b"^] j) and gsib"^] j') in deeply inelastic 
scattering are related to the analogous functions in the other two processes. Consider first 
gi{b), the coefficient of \og{Q'^/Ql). According to the analysis of Ref. [^, this function 
receives contributions from the two jet subdiagrams in Fig. 0(b). (In this figure, we use 
a space-like axial gauge.) The soft gluon connections in Fig. ^(b) affect gA{x,b'^;j) and 
gsip'^]]'), but do not contribute "double logarithms," and thus do not affect giih). Thus 



g^{b)^g?^%b)=gT{b)+gr'\b) 



,OUT/ 



(58) 



where 5'™(&) is associated with the incoming beam jet (the lower subdiagram in Fig. |^(b)) 
while g^^'^{b) is associated with the outgoing struck-quark jet (the upper subdiagram in 
Fig. |^(b)). In the Drell-Yan process, depicted in Fig. |^(a), there are two incoming beam jets 
and one has 

g'^^{b) = 2gf{b) . (59) 
In e^e~ annihilation, depicted in Fig. |^(c), there are two outgoing quark jets and one has 

gf{b) = 2<7?UT(^) _ (60) 

Thus 

g,{b) ^ g^^'ib) = (1/2) gf^ib) + (1/2) gf{b) . (61) 



In the following section, we show numerical results using Ref. ||2^ for gl'^{b) and Ref. 
for gf^{b). 
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0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



b (1/GeV) b (1/GeV) 

FIG. 10. Interpolation of the non-perturbative function e~^^^^^^ vs. 6 as a function of the 
t-parameter, {t = 0, 1/4, 1/2,3/4, 1}, for Q = 30GeV [Fig. (a)], and Q = 100 GeV [Fig. (b)]. Note 
the variation of e~^'^P^^^ as t ranges over [0, 1] is narrower than the full range between the Drell-Yan 
and e+e~ case, (c/.. Sec. [V| ). 



The situation for gA{x,b'^]j) and Qsip'^]]') is not so simple. Let us write 

for the Drell-Yan process and 

^^'p(a;,Q^62;J■ j') = \og(Q' /Qf) gfib) + gf{b'- j) + gf{b'- f) . (63) 



for the energy-energy correlation function in e^e~ annihilation. (C/., Fig. ^.) One might 
like to assume that gA{x,b'^;j) is the same function as g2^ {x,b'^] j) while gsip'^]]') is the 
same function as 5'2^(&^; j')- However, this may not be true because all of these functions get 
contributions from the soft gluon exchanges that link the two jets in Fig. ^ (represented by 
the function U{b) in Ref. |]2S|). Furthermore, the dependence of the functions g^^ {x,h'^] j) 
and gf^ib"^] j') on the flavors j and j' has not been determined from experimental data. What 
we know are flavor averaged functions (7°^(x, 6^) and gl^ip'^)- Thus the best we can do is 
propose a model for the functions we need: 

gA{x,h'-3) + gB{h'-,3') = tg^^{x,b') + il-t)gTib') , (64) 



where < t < 1, with g2^{b'^) taken from Ref. [^ and g2^{x,b'^) taken from Ref. [jl^. We 
vary the parameter t between and 1 to get an estimate of the uncertainty involved. ((7/., 
Fig. m) 



For comparison, we present the above parameterizations for the non-perturbative contri- 
butions with the recent fit by Ladinsky and Yuan for W-production in Fig. |TT]. Ladinsky 
and Yuan introduce an extra degree of freedom by allowing for a, t = xaXb dependence. 
We present the comparison for a range of r; this allows one to gauge the effects of different 
non-perturbative estimates, and correlate the Ladinsky and Yuan parameterization with 
that presented in Eq. (|56D and Eq. 
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b (1/GeV) b (1/GeV) 

FIG. 11. Comparison of the non-perturbative function vs. b for the case of Drell-Yan 

(Davies, Webber, and Stirling |]lj]) [upper dashed Hne], e+e~ (Collins and Soper [^) [lower dashed 
line], and VF-production (Ladinsky and Yuan [^) [5 solid lines]. Fig. (a) is for Q = lOGeV (left), 
and Fig. (b) is for Q = 100 GeV (right). The fits to VF-production include an extra parameter 
T = s/s; we allow r to range over the values r = {10^^, 10^^-^^, 10^^-^, 10^^-^^, lO^^ *^} where 
r = 10~^ is the upper curve, and r = 10"'^ is the lower curve in 6-space. 

VII. REPRISE 

For the benefit of the reader, we review the principal steps in the calculation of the 
energy distribution. The energy distribution is given by: 

dx dQ'^ dqj: d(p ^ dx dQ"^ dqj, dcj) 

= Aki^A) E E MQ';Vi,V2,j,f,k) T,{x,Q\ql-3,f) , (65) 

k=l Vi,V2 

where Ak{ip-, 4>) are the nine angular functions arising from hyperbolic D^{ip, 0) rotation ma- 
trices. The sum on Vi and V2 runs over vector boson types, {7, Z} or {VT^} as appropriate. 
The sums over j and j' include all quark flavors, {u,u,d,d, . . .}; for neutral currents, this 
sum is diagonal (j = j'). The function So(Q^; Vi, V2, j, j', fc) includes factors for the cou- 
pling of the electron to the vector bosons as well as factors for the propagation of the vector 
bosons. The energy distribution function that we have computed is Tk{x,Q'^,q^;j,j'). 

In the limit 0, the Fi and Fg will contain the dominant singularities as their angular 

structure is proportional to the Born process. We define: 

Ti (x, Q^ g^; J, /) = Ff-*(x, Q^ g^; J, /) 

+ r (^1^ { lV(x,Q^g|;J,/)-^f^"^(a;,Q^g^;J,/)} 
^6(x,Q^g^;J,/) = T^''\x,Q\q'T;j,f) 

+ r(^|){(-l)Iy(x,Q^g^;J,/)-Ff^-(a:,Q^g^;J,/)} (66) 
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where the matching function T^qx/Q) [Eq. (|^)] is provided to ensure proper behavior as 
(It Q- rf represents the perturbative results of paper I [Eq. (0)] calculated at order 



a 



1 -pAsj/m j^gpj^gggnts the asymptotic limit {qt — > 0) of Ff'^'"* [Eq. (pi])], and W{x^ Q^, Qt] ]■,]') 



; k 



represents the summed (Sudakov) term [Eq. (p3D ] which is finite as qx — >■ 0. Note the function 
W{x,Q'^,qx]j,j') the same for both Fi and Fg. 

The form of the Sudakov structure function is particularly simple in impact parameter 
space: 



W{x,Q',ql,j,j' 



(2^ 



J, J 



(67) 



To ensure that the calculation is reliable for large b (small qx), we introduce: 

W{x,Q\b';j,j') = W{x,Q',bl;j,j') e-'--^'"^''^^^'^'^ , (68) 

where 6* G [0, bmax] for b E [0, oo]. 

The perturbative function W{x,Q'^,bl]j,j') is given by: 



W{x,Q^bl-j,j') 



Y: fa/A{C,f^) C^ai^Xl^) / dzZ ^ C^.f^Z^ e 



-5(6*) 



(69) 



where x = x/^. For the incoming particles, there is an integration over a parton mo- 
mentum fraction ^, a sum over parton types a = g,u,u,d,d, . . ., a parton distribution 
function fa/ a and a set of perturbative coefficients C™. For the outgoing partons, there is 
an integration over parton momentum fraction z, weighted by z, a sum over parton types 
a' = g, u, u, d,d, . . ., and there are perturbative coefficients C*^^'^ associated with the outgo- 
ing states. The heart of the formula is the Sudakov factor exp[— ^(fo^)], defined as: 



S{b, 





|ln 




Icl/bl /i^ 







A{as{^i)) + B{as{^i)) 



(70) 



The functions A, B, as well as and C^^""", have perturbative expansions in powers of 
ttg. We choose the arbitrary constants {Ci, C2} as in Eq. (^). 

The non-perturbative contribution is parameterized in terms of the fits to e+e" and 
Drell-Yan data.lllilB 



Snp{x,Q ,6 ; j, j ) = log 



Ql 



9?^ib)+gfib) 



+ tg^^ix,l)') + il-t)gl^{i;') . (71) 



The arbitrary parameter t G [0, 1] interpolates between the e~^e and Drell-Yan form. 



VIII. RESULTS 



We present numerical results of the energy distribution function for representative values 
of {x,Q^} using the CTEQ3 parton distributions. We present results only for the Fi set 
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qj (GeV) 

FIG. 12. The contributions to the energy distribution function qj^ dT,i / {dx dQ^ dqj^ d(j)) as a 
function of qt, for Q = 100 GeV, x = 0.3. (RecaU, dSi and dEg are independent of (p.) Perturbative 
(thin), asymptotic (dashed), and summed (thick). Note how the perturbative and asymptotic 
cancel as — > 0. For qr ^ Q, the asymptotic and summed cancel to leading order only. (A zero 
reference line is indicated.) dSi is in units of GeV~^, and is multiplied by 10^ for clarity of the 
plot. 



of structure functions, as the Fg set have the identical qt structure (up to a sign). 
Recall that the structure functions are given by: 

^l(a;,g^g^;J,/) = Ff-*(x,Q^g^;J,/) 

+ ^(f ) {^(^'Q''^T;j,/)-rf'^'"(^,Q',4;j,/)} • (72) 

Making use of Eq. (^, we have a parallel relation for the energy distribution function: 

6;Sl(x,Q^g^;J,/) _ ciSf-^(x, Q^ j,/) 
dx dQ'^ dqj: dip dx dQ"^ dqj< dcf) 

I ^/gT^i c^Sf-(x,Q^g^;J,/) rfsf ^-(o:, Q^ g^; j,/) ] ^^^^ 
\ Q / I dx dQ"^ dq^ dcj) dx dQ"^ dq^ dcp J ' 

where we use the "Sum" superscript to denote the summed Sudakov contribution derived 
from W . We will examine both the individual terms as well as the total in the following. 
We will use the shorthand dSi = dSi (x, Q^, q^] j, j') / {dx dQ^ dq^ dcp) 



A. qT Distributions 

In Fig. |1^ and Fig. |T3|, we show the separate contributions to dT,i as a function of 
for two choices of {x, We have included an extra factor of q^ to make the features 



^ In the small qx region, dSi and dT,Q are independent of (j); therefore we need not specify it. 
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FIG. 13. The contributions to the energy distribution function Qj^ dT,i / (dx dQ'^ dqj^ d(j)) as a 
function of qt, for Q = 30 GeV, x = 0.1. (Recall, dSi and dSg are independent of (p.) Perturbative 
(thin), asymptotic (dashed), and summed (thick). Note how the perturbative and asymptotic 
cancel as qt ^ 0. For Qt ^ Q, the asymptotic and summed cancel to leading order only. (A zero 
reference line is indicated.) dUi is in units of GeV~^, and is multiplied by 10^ for clarity of the 
plot. 



of the plot more legible. As anticipated, we see that dSf ~ rfSf ^™ as — i> leaving 
dSi ^ (iSf""^. For large qt, we find dEf""^ — dsf'^^™ ~ 0, but this cancellation is not 
as precise as the above because the relation F'^"™' — F'^'^^™ ~ holds only to first-order. 
Therefore, in the following figures we shall include the T{q^/Q'^) factor to ensure that 
dYjSum _ dYj^'^y'^ ig smoothly turned off at large qt- The fact that ciSf"™' and c/sf*^"^ 
become negative for large reminds us that these expressions were approximations valid 
only for qt ^ Q- 

Having examined the separate terms, we now turn our attention to the energy distribution 
function, c/Si. Again, we have included an extra factor of in Fig. 0(a) and Fig. p!5|(a) 
to make the features of the plot more legible. In Fig. |l5(b) and Fig. 0(b), we plot dJ^i in 
the small qt region (without an extra factor) to demonstrate that the summed results 
approach a finite limit as gr —>■ 0. We present the results for three choices of the non- 
perturbative function S'np(x, QL,^^; j, j') as parameterized in Eq. (|64D. The choice t = 
corresponds to the e^e~ limit, il while t = 1 corresponds to the Drell-Yan limit, and 
t = 1/2 corresponds to an even mix of the above. The difference due to the non-perturbative 
contribution is quite significant for low qt- The t = {e~^e~) non-perturbative function, 
which is much narrower in 6-space, yields a broader energy distribution; this is clearly 
evident in the figures as we see the peak move to lower qt values as we shift from the 
t = {e~^e~) to t = 1 (Drell-Yan). At large qt, d'Ei is independent of the non-perturbative 
contributions, since it is dominated by dSf 

Clearly, the HERA data should be able to distinguish between this range of distributions, 
particularly in the small qt regime where the span of the non-perturbative contributions are 
significant. 011 
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FIG. 14. The total contribution to the energy distribution function dT^i / {dx dQ'^ dqj^ dcp) as 
a function of qt for different choices of the non-perturbative function, S]\fp{b), for Q = 100 GeV, 
X = 0.3. Fig. (a) has an extra factor of q'^ to make the plot more legible. Fig. (b) demonstrates 
that the summed contribution has a finite limit as qt 0. We vary the t-parameter from t = 1 
(thick) corresponding to the Drell-Yan case, to t = 1/2 (dashed) corresponding to the mixed case, 
to t = (thin) corresponding to the e+e" case. For qx Q, we use the function T^qx/Q) with 
yO = 5 to smoothly switch between large and small qx- dTji is in units of GeV~^, and is multiplied 
by 10^ for clarity of the plot. 

IX. CONCLUSIONS 

Measurement of the distribution of hadronic energy in the final state in deeply inelastic 
electron scattering at HERA can provide a good test of our understanding of perturbative 
QCD. Furthermore, we can probe non-perturbative physics because the the energy distribu- 
tion functions are sensitive to the non-perturbative Sudakov form factor S^p^b) in the small 
Qt region. 

We have evaluated the energy distribution function for finite transverse momentum qx 
at order as in paper I. Because the distribution is weighted by the final state hadron energy, 
this physical observable is infrared safe, and independent of the decay distribution functions. 
In this paper, we sum the soft gluon radiation into a Sudakov form factor to evaluate the 
energy distribution function in the small qx limit. By matching the small and large qx 
regions, we obtain a complete description throughout the kinematic range. This result is 
significant phenomenologically as a the bulk of the events occur at small qx values, where 
perturbation theory by itself is divergent. This technique can provide an incisive tool for the 
study of deeply inelastic scattering. Additionally, crossing relations allow us to relate the 
non-perturbative contribution in deeply inelastic scattering energy distributions to analogous 
quantities in the Drell-Yan and e^e~ annihilation processes. 
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FIG. 15. The total contribution to the energy distribution function dT,i / {dx dQ'^ dq^d(p) as a 
function of qx (scaled by 10^), for different choices of the non-perturbative function, SNp{b), for 
Q = 30 GeV, X = 0.1. Fig. (a) has an extra factor of qj^ to make the plot more legible. Fig. (b) 
demonstrates that the summed contribution has a finite limit as g-r — ^ 0. We vary the t-parameter 
from t = 1 (thick) corresponding to the Drell-Yan case, to t = 1/2 (dashed) corresponding to the 
mixed case, to t = (thin) corresponding to the e^e~ case. For qj- Q, we use the function 
T{qT/Q) with /9 = 5 to smoothly switch between large and small qx- dT,i is in units of GeV~^, 
and is multiplied by 10^ for clarity of the plot. 
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APPENDIX A: KINEMATIC RELATIONS 



We present some basic kinematic relations to facilitate the calculation. First we give the 
expressions to relate {E', 9'} to {x, Q^}, 



= -q.q = 2EE'{l-cos9') 

g2 EE'{l-cose') 



X 



2q-PA Ea[2E - E'il + cosO')] 
Next, we give the expression for the Born scattering angle 6^,, 

1/2 



(Al) 
(A2) 



cot 



\ 2xEa 




1 Q 


xs 



(A3) 



The corresponding azimuthal angle, 0^, is trivial, and can be defined to be zero. Finally, we 
give the expressions to compute the natural variables of the Breit frame, {qx-, 0}: 



. 8^2 _ 4^/ 2^ _ ^' 1 + cos ^' { 9 
qi = ^ < sm^ 

l-cos^B 1 



Ob — d* 



+ sin 6*5 sin6'^. sin^ 



(A4) 
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APPENDIX B: ENERGY DISTRIBUTION FORMULAS 



We now give some explicit formulas for computation of the structure functions and 
energy distribution contributions. The process we consider is the hadronic process e~ + A — > 
e~ + B + X, and the fundamental formula for computation of the structure functions and 
energy distribution contributions is: 



dx dQ'^ dq^ ( 



with 



Y: Aki^PA) E E MQ';Vi,V2,j,f,k)T,{x,Q',q'T;j,j') 

k=l VuV2 



■ (Bl) 



MQ';Vi,V2,j,j',k) 



Q' GUV,,V2;j,f) Gi{V^,V2) 



(B2) 



Ak{ip, 4>) represents the nine angular functions arising from the hyperbolic ^D{iIj, (p) rotation 
matrices. Gl{Vi, V2; j, j') and GliVi, V2) are the combinations of couplings from the leptonic 
and hadronic tensors, respectively, as defined in paper I. (Q^ + My_) arise from the boson 
propagators, and Tk{x,Q'^,q^;j,j') are the hadronic energy distribution function. We sum 
over the intermediate vector bosons {Vi,V2} = {7, Z°} or {IV^}, as appropriate, and the 
parton species 





) = (+1) 




[1 + cosh2(V')] 


A2 (p 


) = (-2) 






A^ {ip, (p 


) = (-1) 


cos(0) 


sinh(2^/') 
sinh^(-?/') 


A (V^, (p 


) = (+1) 


cos(20) 


A^ (V^, (p 


) = (+2) 


sin(0) 


smh.{ip) 


Aq (^, <p 


) = (+2) 




cosh(-0) 


At {ip,(p 


) = (-2) 


cos(0) 


smh.{tp) 


As {ip, (p 


) = (-1) 


sin(0) 


sinh(2?/') 
sinh (-0) 


Aq {ip, (p 


) = (+1) 


sin (20) 



Note, for instance, the analogy between the angular coefficient Ai = 1 + cosh^(?/^), which 
appears in the order energy distribution, and the corresponding coefficient in the case of 
the Drell-Yan energy correlation, 1 + cos {9). [§ 



APPENDIX C: DAVIES, WEBBER, & STIRLING PARAMETRIZATION 



The form of the non-perturbative Sudakov function Spff>{b), used by Davies, Webber, and 
Stirling to introduce the transverse momentum smearing in the Drell-Yan process is: 
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Sfjp{b) 



9i + 92 In 



(CI) 



with 



9i 

92 



0.15 GcV^ 
0.40 GeV^ 



(C2) 
(C3) 
(C4) 



APPENDIX D: COLLINS & SOPER PARAMETRIZATION 



The form of the non-perturbative function used by CoUins and Soper to introduce the 
transverse momentum smearing in the e"*" e~ process is: 



5^(6)= A {4^1^ In 



TT 



C2QK 



m - +AM5)ln 



+ AMb) (Dl) 



with 



AMb) = Anb + Aub^ 
A/2(6) =^21 6 + ^22 



(D2) 



While the functional form allowed here is quite general, in practice, it was possible to obtain 
a good fit to the data using only the A and A21 parameters. Specifically, 



A = 1.33 
A21 = 1.5 

An = Au = ^2 = 

Additional parameters and relations necessary are: 

&ma. = (2GeV)-^ 
Qo = 27GeV 

Ai^2Cf 



(D3) 



(D4) 



APPENDIX E: LADINSKY & YUAN PARAMETRIZATION 

The form of the non-perturbative Sudakov function S^pib), used by Ladinsky and Yuan 
to introduce the transverse momentum smearing in the Drell-Yan process is: 



SNp{b) 



9ib'^ + 91 93bln[100T] + g2b^ In 



2go. 



(El) 
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with 



gi = 0.11 GeV^ 
g2 = 0.58 GeV^ 

Qo = 1.60 GeV^ 
6^„, = (2GeV)-i . (E2) 



APPENDIX F: as AT 1-LOOP AND 2-LOOP 

To properly compute the fi^ integral in the Sudakov form factor, it will be necessary to 
use the complete result for the running coupling at both 1- and 2-loops. The 2-loop result 
for a, is: 



where 



_ 47r 47r/32ln[ln(/xVA=^)j 



^ inK-2Nf) ^ (33 - 2Nf) 
^^3 3 ^ ' 

ft = (102 -5^) . (F3) 



The 1-loop result is simply obtained by taking (32 —>■ 0. 



APPENDIX G: INTEGRAL TABLE 

For simphcity and completeness, we list the integrals we shall encounter in the Sudakov 
form factor at the 1- and 2-loop level. We consider the logarithmic terms (Ai) and the 
constant terms (Bi) using the 2-loop expression for a^; the 1-loop expressions are easily re 
covered in the limit /?2 0. It will be convenient to define the following quantities: 



Li = In 



Lo = In 



In 



62 A2 

CI 

A2 



= Li — L3 



(Gl) 



First, the Ai term with the 2-loop expression for ct^ 
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In 



IJ? \ (2)7r 



L2 + L3 In 



Li 



+ 



The Bi term with the 2- loop expression for a^. 



Bi — „ In 



clihl tJ? (2)7r 



+ 



(2)/3i 
4/^25 



^1 



(2)/33LiL2 

The term with the 1-loop expression for ag- 



i-(^i-i^3 + i^iln[L3]-L3ln[Li]) 



L 



clQ^ df^ (ajji]iy 
cf/bi IJ? \ (2)7r ^ 



A, 



16A 



— L2 — Li In 



^1 



APPENDIX H: BOSON-FERMION COUPLINGS 



Fermions 




9a{l) 


9viZ) 











— e 


1-4 sin^ 6yy 


+e 


1 


— e 


4 cos Spv? sin Ow 


4 cos 9w SID 0w 


tt, c, t 







+e 


l-|sin2 6»vy 


— e 


1 




4 cos Ow sin ^v^/ 


4 cos sin 




-I- 





— e 


1 — 1 sin^ 9w 

12 


+e 


1 


4 cos Ow sin Svv' 


4 cos Svi/ sin 6w 



Table 1. Boson-fermion couplings. 
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